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1 Introduction 

Discrete symmetries play a fundamental role in particle physics. Charge conjugation (C) 
and the discrete spacetime symmetries of parity {V) and time reversal (T) are such that 
CVT is necessarily conserved for a local, Lorentz-symmetric, and Hermitian theory. There 
is, however, no reason for all of the latter requirements to be essential in the building of 
viable models. This has been shown, for example, in ref. [1], where locality is dropped, 
leading to a Lorentz-symmetric description of neutrino physics, which is odd under CVT. 
In the present article we keep locality but consider a non-Hermitian (Lorentz-symmetric) 
model. 

The last 15 years have seen much interest in and research activity on theories described 
by non-Hermitian Hamiltonians. Such theories have remarkable and often unexpected 
properties. For example, the eigenvalues of the non-Hermitian CPT-symmetric quantum- 
mechanical Hamiltonians H = + ix^ and H = p^ — are real, positive, and discrete 

[2, 3]. 

One idea that has been pursued repeatedly is to study the properties of a non-Hermitian 
version of quantum electrodynamics (QED). The Hamiltonian for QED becomes non- 
Hermitian if the unrenormalized electric charge e is chosen to be imaginary. Then, if 
the electric potential is chosen to transform as a pseudovector rather than a vector, the 
Hamiltonian becomes CVT symmetric. The resulting non-Hermitian theory of electro¬ 
dynamics becomes a multi-component analog of a self-interacting spinless quantum field 
theory (QFT), comprising a pseudoscalar field (j) with a cubic self-interaction term of the 
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form . This pseudoscalar QFT was studied in detail in ref. [4], and the non-Hermitian 
version of QED was studied in ref. [5]. This non-Hermitian version of electrodynamics is 
particularly interesting because it is asymptotically free and the version of this theory with 
massless fermions appears to have a nontrivial fixed point (see refs. [6-8]). A perturbative 
calculation of a metric with respect to which this theory is unitary is given in ref. [9]. 

A detailed analysis of CPT-symmetric non-Hermitian fermionic theories was done by 
Jones-Smith and Mathur [10]. In this work it was emphasized that for fermions the time- 
reversal operator T has the property that = — 1. This represents a significant departure 
from the case of bosonic theories, where = 1. (Further work on the properties of CVT- 
symmetric representations of fermionic algebras may be found in ref. [11].) In addition, 
Jones-Smith and Mathur showed that free noninteracting CPT-symmetric Dirac equations 
have the remarkable feature that massless neutrinos can exhibit species oscillations [12]. 

The discovery of neutrino oscillations and the observation of the baryon asymmetry 
of the universe (BAU) (see ref. [13]) have been driving forces in the study of the neutrino 
sector of the SM. Neutrino oscillations consistent with experimental observations can occur 
if the SM neutrinos have small but finite masses. The misalignment of the mass and 
flavour eigenbases then gives rise to the PMNS [14, 15] mixing matrix, analogous to the 
CKM [16, 17] mixing matrix of the quark sector. In order to generate the BAU, it is 
necessary to satisfy the Sakharov conditions [18]: namely the presence of out-of-equilibrium 
dynamics and the violation of baryon number B, charge C, and charge-parity CV. Both 
the CKM and PMNS matrices contain a complex phase, which provides a source of CV 
violation in the SM. In the quark sector this gives rise to the CV violation observed in 
K-, D-, B- and H^-meson mixing (see ref. [13]). However, the magnitude of this CV 
violation is insufficient to have generated the observed BAU. An elegant framework in 
which both experimental observations may be accommodated is provided by the scenario 
of leptogenesis [19] (for reviews, see refs. [20-23]). Therein, the SM is supplemented with 
heavy Majorana neutrinos. The smallness of the light neutrino masses arises by means of 
the see-saw mechanism [24-28] and the baryon asymmetry through the decays of the heavy 
neutrinos in the expanding early universe. By virtue of the lepton-number L-violating 
Majorana mass terms and complex Yukawa couplings, which provide an additional source 
of C and CV violation, these decays are able to generate an initial lepton excess, which is 
subsequently converted to a baryon excess via the {B -|- L)-violating electroweak-sphaleron 
interactions of the SM [29]. 

In this article we examine an extension of QED that involves the usual Dirac mass 
term and an anti-Hermitian mass term The fermion field is coupled to the 

photon through both vector and axial-vector couplings. The anti-Hermitian mass term is 
separately C even, V odd and T even, and is consistent with unitarity for /i^ < m^. We 
study the gauge symmetry of this model and show that, although gauge invariance is lost 
in the massive case, it is recovered in the specific situation where the Hermitian and anti- 
Hermitian mass terms have equal amplitude = m^. In this limit we find that the model 
reduces to that of a massless left- or right-chiral Weyl fermion. Moreover, we illustrate 
that by choosing the ratio n/m we may obtain an arbitrarily small but finite mass for the 
fermion and give more or less prominence to one chirality. This observation, combined 
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with the maximal CV violation of the anti-Hermitian mass term, may be directly relevant 
to neutrino physics. 

The paper is organized as follows: Section 2 begins by summarizing the essential 
properties of the free non-Hermitian fermion theory studied already in refs. [30] and [31]. 
Subsequently, the gauge interactions are introduced and the tree-level properties of the 
model are described. Therein, emphasis is given to the restoration of gauge invariance 
in the limit jj? = m^. Section 3 presents the one-loop self-energy and vertex corrections, 
the details of which are given in appendix A. Here, the recovery of gauge invariance is 
made explicit through the expected vanishing of the longitudinal component of the vacuum 
polarization. Section 4 describes an analogous non-Hermitian Yukawa model and discusses 
possible implications for the neutrino sector of the SM. A novel mechanism for generating 
the light neutrino masses as well as the presence of an additional source of CV violation is 
highlighted. Concluding remarks are given in section 5. 

2 Description of the Model 
2.1 General description 

We begin with the free fermion non-Hermitian Lagrangian considered in ref. [30]: 

Cq = {i$ - m - /i 7 ®) il) , (2.1) 

with fj? < m^, such that the energies uj are real for all three-momenta p; that is, 

a ;2 = ^2 > 0 , ( 2 . 2 ) 

where 

_ ^2.3) 

It is shown in ref. [31] that the conserved current for this model is 

jp = iPY U + (2.4) 

V m / 

and that the equation of motion is obtained by taking the variation of the action with re¬ 
spect to Ip for fixed ip. The anti-Hermitian mass term in eq. (2.1) is even under both charge 
conjugation C and time-reversal T, and odd under parity V. Thus, it is odd under CVT. 
However, this does not contradict invariance under Lorentz-symmetry, since Hermiticity 
has been relaxed. 

In this article we gauge this model and include both vector and axial-vector coupling 
to an Abelian U{1) gauge field A^: 

£ = - I + il)[i$ - ^{gv + gAl^) - m - ip , (2.5) 

where F^y = d^Ay — dyA^. In the massless case m = g = 0 the action is invariant under 
the combined vector and axial gauge transformation 


\ - 

A^ d^(p , 

(2.6a) 

Ip - 

exp [i [gv + gAl^) 0] V’ , 

(2.6b) 

Ip - 

Ip exp [i (- gv + gAl^) (p] ■ 

(2.6c) 
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However, in the massive case m 7 ^ 0 and/or fi ^ 0 this gauge invariance is lost. 
The free fermion propagator of this theory is 

, # + m - ^ 7 ® 

*‘3 = * -w , 

+ l£ 


(2.7) 


where e = O"*". We see immediately that eq. (2.7) has a light-like pole for /r = ± m (M^ = 0), 
like that of a massless theory, with the propagator taking the form 


iS = 



( 2 . 8 ) 


The mass term in the numerator is proportional to the chiral projection operators 




(2.9) 


where is the n x n unit matrix. Separating the right- and left-chiral components tpR = 
Pr'iP and 'tpi = Piip in the current (2.4), we see that the probability density may be written 
as 

p = (1 + -) IV'rI" + (1 - IV’Lp . (2.10) 

V mJ V mJ 

Evidently, for p = +{—)m the contribution to the probability density is entirely from the 
right-(left)-handed degree of freedom. Therefore, it appears that in the limit /x = -|-m we 
obtain a massless right-handed theory, and in the limit p = —m we obtain a massless left- 
handed theory. This feature is the focus of this article. Moreover, in Sec. 4, we comment 
on potential implications of this non-Hermitian theory for the neutrino sector of the SM 
and, in particular, the smallness of the light-neutrino masses. 

The preceding observations suggest that it proves illustrative to consider this theory 
in an explicit chiral basis. We do so in the following section and show explicitly that 
invariance under the gauge transformation in eq. ( 2 . 6 ) is recovered in the limit p ^ Pm, 
as we would anticipate for a theory that appears to be effectively massless. 


2.2 Chiral basis 

In order to recast eq. (2.5) in an explicit chiral basis, we first rotate from the Dirac basis 
to the Weyl basis via the orthogonal transformation 


'IpW 




( 2 . 11 ) 


We may then work directly with the two-component right- and left-chiral spinors ipR and 
IpL- 

In the Weyl basis the gamma matrices take the form 


n 

Iw — 


0 

0 


Iw = 


(-I2 o\ 

\ 0 hj 


( 2 . 12 ) 
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where = ((T°,(T®) and = (cr'^, — it*), and u* are the Pauli matrices. To avoid a pro¬ 
liferation of subscripts and superscripts, throughout this paper we suppress SU ( 2 ) spinor 
indices (see appendix A). In addition, the projection operators are given by 


p. = r*" 

The fermionic sector Lagrangian is then 
•^ferm = V’Jj 

where 


with 



)■ 

(2.13) 

/ ia ■ D- — \ / 

1 — m_ ia ■ D_|_ j \ 

't) ■ 

,’PrJ 

(2.14) 

= m ± ^ , 


(2.15) 

+ ig^A^ , 


(2.16) 

= 9V ^ 9A ■ 


(2.17) 


Notice that 7 ® matrices nolonger appear explicitly in the Lagrangian eq. (2.14). Instead, 
the non-Hermitian nature of this theory is manifest in the asymmetry between the right- 
and left-chiral components of the original four-component Dirac spinor. 

We may study the on-shell structure of the Lagrangian in eq. (2.14). For the case 
fj, = +m the Lagrangian takes the form 


= + ^ = 'iplia ■ D.-ipL + ■ D+'ijjR - , (2-18) 

giving the following equations of motion for 'i/jr and 'ipi: 

s s 

— T- = 0 => ia ■ Dj^TpR = 0 , (2.19a) 

dV-L 
5 s 

—Y = 0 ia ■ D-TpR = 2m'4)R . (2.19b) 

Since the left-chiral field does not appear in the equation of motion for the right-chiral field 

[eq. (2.19a)], we may integrate it out, giving the tree-level on-shell Lagrangian 

Cn-shell = ■ D+^JR , (2.20) 


which describes a massless theory of right-handed Weyl fermions. This is precisely what we 
saw in the probability density [eq. (2.10)]. Moreover, the on-shell Lagrangian [eq. (2.20)] 
respects the full vector and axial-vector gauge invariance [see eq. ( 2 . 6 )]; that is. 


An - dn 




exp (ifi(+(/>) 'il)R 


( 2 . 21 ) 


For the case /r = —m we need only make the replacements 'i/jr -h- 'ipi, d -f-)- a, and 
D+ -f-)- D_ in eq. (2.20), yielding a massless theory of left-handed Weyl fermions. The next 
subsection gives a more explicit argument to justify the restoration of gauge invariance for 
the light-like case fj? = m?. 
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2.3 Hidden gauge invariance 

In this subsection we show that gauge invariance is recovered when fj? = To do so, we 
construct a two-component spinor basis in which gauge invariance is explicit. 

Written in block form, where the LL (left-left) element is located in the top left 2x2 
block, the mass matrix is given by 


M = 


0 m+ 
m_ 0 


( 2 . 22 ) 


having eigenvalues ±M = — iJ? and eigendirections 


e± = 


^ with x± = ± /r/m . 


(2.23) 


We rotate to the mass eigenbasis but first allow for a rescaling of the left- and right-handed 
components: 

V'L(R) —> V'l(r) = aL(R)'0L(R) , (2.24) 

where to be determined later, as explained below. This leads to the transforma¬ 

tion 


V’+) _ .^-1 ( V’L 


>- 




with 


v- = T 


aLX- aRX+ 

V2 y-aix- aRX+ 

The Lagrangian then becomes 


and V = T ( 

V2 \ l/(aKX+) l/{aRX+) 


>Cferm = V’l) ^ 

= V’l 


ia • Z)_ — m+ 
— m_ ia • 


V 


i’+ 


A B 
C D 


'4^+ 


where 


ia • Dj^ 

+ 

\ 

lb 


m 

2a%xl 

2a|x?_ 


aRttLX+X- 

ia • 


ia • D- 



2alxl 


2a\x‘^_ 


aRttLX+X- 

ia ■ Dj^ 


ia • 

+ 


2alxl 


2alxl 

aRttLX+X- 

ia ■ Dj^ 

+ 

ia • D- 

+ 

m 

2a\xl 

2a\x‘^_ 

aRttLX+X- 


(2.25) 


(2.26) 


(2.27) 

(2.28a) 

(2.28b) 

(2.28c) 

(2.28d) 


- 6 - 



and the two-component spinors V’-i- and '!/>_ are given by 

V'i = ± . (2.29) 

The next step is to determine the coefficients (1l(r)- To make gauge invariance explicit 
in the limit /r —)■ +(—)m, that is, —)■ 0, only the covariant derivative should 

remain. A reasonable choice for the field rescaling is 


X' 


2a|x2 


= ^:r and 


2a\x‘^_ 


9 ’ 


where a is an overall numerical coefficient. Thus, we have 

a , a 


CLR = 


V2xl 


and 


ai = 


V2xl ’ 


(2.30) 


(2.31) 


and we obtain 


^ -^ferm — 


1 'j/d 7/d 1 

1 x\ia ■ 

• + x^_i(J ■ 

• T»_ 

x\ia ■ 

■ D+ — x‘^_ia ■ 



1 ' 

• — x^_ia ■ 

• T»_ 

x\ia ■ 

■ Dj^ -|- x‘^ia ■ 

■ d) 





(2.32) 


Note that the mass matrix is not diagonal, even in the mass eigenbasis, because of the 
anti-Hermitian mass term controlled by /i. 

In the limit /U —)• ± m the mass term vanishes, and we are left with a massless theory 
that is invariant under the gauge transformation 


An - dn4> 


V'li 


exp (ic/+(/>) V'i , 
exp (ifif_</>) V'i , 


Moreover, from eq. (2.29), we have 


/i = +m, 
/i = — m . 


(2.33) 


iP'r, = +m, 

±V’l , n = -m , 


(2.34) 


and 


-^ferm 


i;lia ■ , 


^ = +m, 
/i = — m , 


(2.35) 


for massless right- and left-chiral theories, as observed in the preceding subsections. 

The coefficients and the transformation in eq. (2.26) are singular in the limit 

// —>• ± m. However, the coefficients do not appear in the final Lagrangian [eq. (2.32)], 
which remains finite in the limit /r —)• ±m. Furthermore, the functional Jacobian of the 
field rescaling in eq. (2.24), although also singular, cancels out in the normalization of Z 
with the partition function Zq of the corresponding free theory. Thus, the limit ^ ^ zizm 
may be taken safely, as done above. 
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2.4 Exceptional points 


An A^-dimensional Hermitian matrix always has N real eigenvalues and associated with 
each eigenvalue is a distinct eigenvector. For non-Hermitian matrices the situation is more 
elaborate. Consider, for example, the non-Hermitian 2-dimensional matrix 


A 


a + ib g 
, g a-ib. 


(2.36) 


where a and b are real parameters and (7 is a real coupling constant. The eigenvalues of 
A are E{g) = a ± \^g^ — b‘^. Thus, there are two phases: a broken phase (for g'^ < 6 ^) 
in which the eigenvalues are complex and an unbroken phase (for g'^ > 6 ^) in which the 
eigenvalues are real. At the boundary between the phases {g = Eh) the eigenvalues merge, 
and the matrix is said to be defective because there is only one eigenvalue E = a and one 
eigenvector instead of two: {i, 1) fox g = + b and (l,i) for = — b. The point g'^ = b^ is 
called an exceptional point. 

In general, at an exceptional point a pair of eigenvalues of a non-Hermitian matrix 
merge, and one of the eigenvectors disappears. (It is possible for more than two eigenvalues 
to merge at an exceptional point, but this is not common.) If the exceptional point occurs 
when a parameter g has the value go, the eigenvalues E{g) exhibit a square-root singularity 
at = go- 


Hermitian matrices do not have exceptional points. Nevertheless, exceptional points 
play a crucial role in explaining their behavior. For example, in conventional Hermitian 
quantum theory the radius of convergence of a perturbation expansion in powers of a 
coupling constant is precisely the distance to the nearest exceptional point (a square-root 
singularity) in the complex-coupling-constant plane [32]. 

In the limit g, = Em the mass matrix M is defective and, as explained above, the 
transformation in eq. (2.26) becomes singular. For instance, for g = Em the mass matrix 
has the Jordan normal form 


M 


2 m 



(2.37) 


with zero eigenvalues. In this limit we have chiral flips biased from left to right, depleting 
the probability density of the left-handed component, as we saw in eq. (2.10). In other 
words, we again arrive at a massless theory of right-handed fields. Conversely, in the limit 
// —)• — m we arrive at a massless theory dominated by left-handed helds with chiral flips 
biased from right to left. 

The appearance of defective matrices is rare in physics, especially in field theory, and 
it is worth considering what this singular behaviour signals. For g'^ < mf we have right- 
and left-chiral components with positive mass-squared > 0 (time-like) and real-valued 
energies a; G M. For g > the right-chiral component dominates; for g < the left-chiral 
component dominates; for = 0 we have exact symmetry between both components. On 
the other hand, for g^ > m? we still have right- and left-chiral components, but these are 
now tachyonic, having negative mass-squared < 0 (space-like) and imaginary-valued 
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Figure 1. Schematic representation of the m-/i plane, where the tachyonic region (cross-hatched) 
is bounded by the lines = m^, along which the mass matrix becomes defective: fj, = + m (red) 
corresponds to total left-chiral domination, and fj, = — m (blue) corresponds to total right-chiral 
domination. Along the line /i = 0 the symmetry between the right- and left-chiral components is 
restored. 

energies iuj £ M. For the special case iJp = m? we have a massless fermion = 0 (light¬ 
like) and real-valued energies lo = \p\ £ M. For p = + m this held is completely dominated 
by its right-chiral component and for /r = — m it is completely dominated by its left-chiral 
component. The mass matrix becomes defective at the boundary between the time-like 
particle and space-like tachyonic regimes. This is indicated graphically in hgure 1. 

3 One-loop corrections 

We give here the one-loop corrections to the fermion and photon self-energies, as well 
as the three-point vertex. The technical details of the calculations are given explicitly in 
appendix A. For our purposes it is convenient to express the one-loop results in terms of the 
Passarino-Veltman form factors [33], the dehnitions of which are also given in appendix A. 
We work in the Feynman gauge throughout, with the gauge-fixing term 

^ : (3-1) 

in which the photon propagator has the simple form 

iD,,{k) = . (3.2) 

+ ie 
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Fermion self-energy. There are four contributions to the one-loop fermion self-energy: 
one with two vector couplings, one with two axial couplings, and two with one vector 
and one axial coupling. Employing dimensional regularization and working in d = 4 — 2e 
dimensions, we find the total self-energy 

^{p) = iav - 9Al^f fBi + -5a)(”^ + ^7'’)-Bo , (3.3) 

where we have suppressed the arguments on the form factors i?o;i = M, 0). Isolating 

the logarithmically-divergent part, we obtain 

^{p) = {^{av + 9A'y^ f+4.{glr - g\){m ++ finite , (3.4) 

where higher orders in p are omitted. 

The RL and LR components of the fermion self-energy are given by 

^Rl{p) = (3.5a) 

Slr(p) = y^dm+Ro- (3.5b) 

For /i = ± m we see that one of these components vanishes such that it remains the case 

that only the operator {p = +m) or {p = —m) is present, preserving the 

argument in Subsec. 2.2. Specifically, the equations of motion for the right- and left-chiral 
fields at order are given by 

• D+iI^r = 0 , (3.6a) 

= (2m + 5m - T.lr{p))' tpR , (3.6b) 


where 5m is the mass counterterm and in the on-shell scheme the wavefunction renormal¬ 
ization is given by 


^R{L) — 1 + ^^RR(LL){P) 


(3.7) 


p^ = o 

We may again integrate out the left-chiral component, obtaining a massless right-handed 
theory also at order g^. 


Polarization tensor. There are also four contributions to the one-loop photon polariza¬ 
tion tensor, and we find the total polarization tensor 

2 2 2 
n'^^(p) = - {P^P" - (^21 + 5i) + ||r/^^M2Ro . (3.8) 

The form factors are evaluated as Ro;i; 2 i = Bo-i- 2 i{p, M, M). As expected from the loss of 
gauge invariance in the case of axially-coupled massive fermions, the polarization tensor is 
not transverse and contains the longitudinal part 

n^^(p) D (3.9) 

TT^ 
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Nevertheless, this longitudinal part vanishes and the polarization tensor becomes transverse 
when fj? = m? (that is, when M = 0). Hence, as a consequence of the restoration of 
gauge invariance (see Subsec. 2.3), the polarization tensor satisfies the standard QED 
Ward identity. 

Isolating the logarithmic divergences in eq. (3.8), we obtain 

2 2 2 

n^-(p) = (pV - + finite , (3.10) 

where higher orders in p are omitted. 


Vertex. The four different contributions to the three-point vertex lead to the total one- 
loop correction 


a ^( p , q ) = \^^3v + gAi^) 


X {(w + [(2 - d)rC2A + + {gl + ffi)7^M2Co} 


1 


+ ^ (9v+5i)(m-p7®) 

X [gv [p^(2Cn + Co) + q^{2C^2 + ^o)] - gAl\p>^ + q^Co} , (3.11) 

where we have defined 

Fnp = PnPp {C\l+C2l) + qK.qp {C22 ~\~Ci2^ + P^qp {C2‘i + C\\^ + q^Pp {C20~\~Ci2) ■ (3.12) 


The three-point form factors (see appendix A) are evaluated with arguments pi = P-,P 2 = q, 
mi = m 3 = M, and m 2 = 0. The divergent contribution to eq. (3.11) arises from the form 
factor C 24 and is given by 

A'" D + 


which is consistent with the self-energy [eq. (3.4)], as imposed by the Ward identity for 
—)■ 0, describing the usual gauge invariance of vectorially-coupled massive QED. 
Einally, the RL and LR components of the vertex, that is, those mediating right-to-left 
and left-to-right chiral flips, are given by 


_ 9 + 9 - 

- 4^2 

_ 9+9- 

- 4^2 


{ 9 + + 9-){p^ Cii + q^ C 12 ) + 9-{p^ + q^)Co] , 
( 9 + + 9-){p^ Cii + q>^ C 12 ) + 9+{p^ + q^)Co 


(3.14a) 

(3.14b) 


Like the RL and LR components of the fermion self-energies, these terms are proportional 
to m_ and m+, respectively, so we have only left-to-right chiral flips for p = +m and 
right-to-left chiral flips for p = —m, which preserves the structure observed in Subsec. 2.2 
for p"^ = w? also at order g^. 
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4 Implications for neutrino masses 

This section highlights potential implications of the behavior of this non-Hermitian theory 
for the neutrino sector of the SM. We extend the SM with a right-handed singlet neutrino 
vr. In the Dirac basis we write the non-Hermitian neutrino Yukawa sector, assuming only 
a single generation for now, as 

C = LlUPLl + VRi$i^R - h-LicpRR - h+VR(j)^LL , (4.1) 

where Lr = is the SU{2) lepton doublet, (j) = ia 2 (j)* is the isospin conjugate of the 

Higgs doublet and is the usual covariant derivative of the SM gauge groups. Note that 
we have swapped -|- and — relative to the non-Hermitian model of QED in the preceding 
sections. The non-Hermitian Yukawa couplings are 


h± = h ± T] , 


(4.2) 


where, for now, we assume that /i, r/ G M. Since the electroweak sector of the SM is already 
written in terms of chiral fields, no 7 ® appears explicitly in the non-Hermitian Lagrangian. 
Even so, in the symmetry-broken phase, the non-Hermitian Yukawa couplings give rise to 
a Hermitian mass m = vh and an anti-Hermitian mass /x = vrj, with m± = v{h ziz rj), in 
complete analogy to the non-Hermitian Abelian theory considered in Sec. 2. 

In the unitary gauge and after spontaneous symmetry breaking, the Higgs doublet 
takes the form 


<P = 


1 

71 





(4.3) 


Hence, the neutrino sector becomes 


C, 


1 ) 

VLi$i'L + VRi$VR - h^^ViVR - h+ 


V - - TT 

7! 7! 


h+ 

%/2 


urHi/l ■ 


(4.4) 


The first four terms of the Lagrangian in eq. (4.4) can be written in the matrix form 




I'L I'R 


where the neutrino mass matrix 



M 


V I 0 h- 

V 2 \h+ 0 


has eigenvalues 

± M = ± V/l2 - 7/2 . 

V 2 

Proceeding in analogy to Subsec. 2.3, we make the field redefinition (a = \/2) 

^ , _ ^L{R) 

^L{R) ^ 1^l{R) - „2 ’ 

^-(+) 


(4.5) 


(4.6) 


(4.7) 


(4.8) 
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where 


(4.9) 


x± = \/l± r]/h . 

We then move to the mass eigenbasis spanned by the two-component spinors 


with 



1 _ 1 


' ^/2 


1 

R ± 

■ 



{x+v'l ± 


(4.10) 


(4.11) 


Thus, in the limit rj ^ h and in analogy to Subsec. 2.3, we obtain a theory of massless 
left-handed neutrinos, which is the “original” Standard Model. However, arranging for 
r/ ~ /i with r] < 0, we obtain a nonzero but arbitrarily small mass for the neutrinos, with 
the propagating state still dominated by its left-chiral component. 

In the above minimal extension of the SM the singlet neutrino ur does not couple to the 
SU{2)l gauge fields, and we cannot make use of an analogy to the non-Hermitian Abelian 
gauge couplings of Subsec. 2.1. However, in the so-called left-right SM [34-36], where the 
SM gauge groups are extended from SU{2)l ® U{1)y to SU{2)l ® SU{2)r ® U{1)b-l, 
the SU{2)r gauge fields couple directly to the right-handed neutrino current. Thus, by 
introducing couplings g+ = gv+SA and g- = gy —9A of the left- and right-handed currents 
to the charged gauge fields and and g'j^ = g'y + g'j^ and g'- = g'y — g'j^ to the 

neutral gauge fields and Z^, those of the right-handed neutrino may be suppressed for 
gry(') ^ gA^'^- This, of course, amounts only to choosing different values for the tree-level 
SU{2)l and SU{2)r gauge couplings, which need not result from a non-Hermitian theory. 
Nevertheless, this construction might provide a common origin for such a structure in the 
gauge and Yukawa sectors. 

The masses of the left- and right-handed neutrinos are degenerate in this construction, 
both being light for rj ^ h. However, since the right-handed neutrino is still a singlet of 
the SM gauge groups, we are not precluded from adding a Majorana mass term 




(4.12) 


where C denotes charge conjugation. In this case, the Lagrangian takes the form 


- Cy ^ 


1 

2 


VL 


( 0 m_ 
m+ niR 




{^L ^r) 


( 0 m_ 
m+ niR 



(4.13) 


For rriR ^ 2M the masses of the light and heavy neutrinos are ttil = — M'^/niR and ttir, 
which drives up the mass of the right-handed neutrino and further suppresses that of the 
left-handed neutrino by means of the see-saw mechanism [24-28]. 

It is worth commenting on the generalization to complex Hermitian and anti-Hermitian 
Yukawa couplings h and g. In this case, the Yukawa sector takes the form 


Cy D 


-h- 


V 

72 




h* 


V 

7^ 


rrrl ■ 


(4.14) 


- 13 - 



The mass-squared is then given by 


2 

“ 2ilm.h*r^ , (4.15) 

which delivers real masses only when Im/i*r/ vanishes, i.e. when h = rj. Thus, if we want 
small but finite masses, we are required to take h and rj to be real. 

The situation is somewhat different, however, when we consider the extension of the 
above model to include N generations: 

C = I’lUpLi^k + i'Ri0i'R,a - , (4.16) 


where we have assumed only Dirac masses in the first instance. We have employed the 
flavour-covariant notation of ref. [37], where the left- and right-handed sectors transform 
in the fundamental representation of two flavour groups Ul{N) and Ur{N), respectively, 
and flavour indices are raised and lowered by complex conjugation. We have taken the 
number of left- and right-handed fields to be equal for simplicity in what follows; this need 
not be the case in general. Under a general transformation in Ul{N) x Ur{N), we have 

LL,k L'^,k = Li = (4.17a) 

l^R,a l^R,a = , I/r = {j^R,a)^ = U^p^^R , (4.17b) 

where V\ = (14 0* ^ Ul{N) and = {Ua^)* G Ur{N). The Yukawa coupling matrices 
h± = /i±r7 transform as tensors of Ul{N) x Ur{N) and flavour covariance of the Lagrangian 
requires the transformation property 

(4.18) 


In general, there will not exist a flavour basis in which the Yukawa matrices and /i_ are 
simultaneously diagonal. As a result, there can be a three-fold misalignment for general 
Yukawa matrices, i.e. the weak, -|- Yukawa and — Yukawa bases can point in three different 
directions in flavour space. Hence, for three generations, neutrino oscillations in this model 
are governed by 6 rather than 3 mixing angles and 2 rather than 1 CP-violating phases. 
This additional source of CV violation is of particular relevance to the potential embedding 
of this non-Hermitian theory within the scenario of leptogenesis. 

In the symmetry-broken phase, for the case of two generations {N = 2), the mass 
spectrum contains four mass eigenstates with masses given by the roots of 


.,2 r 


M4o^ = — 


' 1 ( 2 ) 


4 L 


tr h\h 


— (-I-) ^2tr (h,]^/i_)^ — (tr 


1/21 


(4.19) 


It is clear that one may obtain the massless limit by choosing = ± r/. However, such a 
constraint is not a necessary condition for obtaining a spectrum with massless states. In 
the case that 


deth\_h^ =0 =► tr(/i^_/i_)^ = (ti , (4.20) 
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we obtain two massless states (± Mi = 0 ) and two states with masses given by the roots 
of 

„.2 „,2 

tih^h — — 2iImtrh,^T7 


m 2 = ^tr/rU_ = - 


For M 2 to be real, we require 

Imtr/i^ry = 0 . 

Subsequently imposing the additional constraint that 

tih^h = tr 77^77 , 


(4.21) 

(4.22) 

(4.23) 


we also obtain M 2 = 0, giving a massless spectrum. In complete analogy to the single¬ 
flavour case, we may obtain an arbitrarily small but finite mass splitting AM^ = M| — M^ 
by choosing 

tih^h ~ tr 77^77 . (4.24) 

As a result, there is the potential to obtain sub-eV scale Dirac neutrino masses from Her- 
mitian and anti-Hermitian Yukawa couplings, whose orders of magnitude may themselves 
be much closer to the other SM Yukawa couplings and larger than the unnatural 10 “ 
that would otherwise be required for agreement with neutrino oscillation data. 

Note that eqs. (4.20), (4.22) and (4.23) comprise three constraints on the total of 16 
parameters in the complex-valued 2x2 matrices h and 77 . These three necessary conditions 
provide a much weaker constraint on the elements of h and 77 than the condition h = ij. 
Moreover, they do not, as in the single-flavour case, require h and 77 to be real-valued 
matrices. 

As for the single-flavour case, we can include a Majorana mass term of the form 

1 


- Cy ^ 


2 + H.c. , 

where the mass matrix mu transforms as a rank-2 tensor of Ur{N), i.e. 


a/3 


al3 


In block form the mass terms are given by 
-Cy 


7(5 

R 


(4.25) 


(4.26) 



The mass matrix 


M = 


0 7n_ 

ml mn^ 


(4.27) 


(4.28) 


can be block diagonalized by a unitary transformation of the form M = MW, giving 

the physical neutrinos 


Nl 


= w' 




R/ 



= 



(4.29) 
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where the Nl are the light neutrinos, whose mass matrix is given by the non-Hermitian 
see-saw formula 


rriL = — 

and the Npi are the heavy Majorana neutrinos, whose mass matrix is rriR. 
For N = 2 the mass spectrum of the light neutrinos is given by 


,,2 r 


^1(2) - “ 4 L 


1 / 2 ' 

“ (+) ^2tr — (tr 

We trivially obtain a massless spectrum for h = zizrj. However, as before, when 
deth^mj^^h^ = 0 => tr , 

we obtain the spectrum 


(4.30) 


(4.31) 


(4.32) 


Ml = 0, M 2 


V _ _ 1 , T 

— tr n, 

2 H 


(4.33) 


For M 2 to be real, we now require 

Imtr = 0, 

and we obtain a completely massless spectrum if, in addition, we require that 

Retr= 0 . 


(4.34) 


(4.35) 


Again, the conditions eqs. (4.32), (4.34) and (4.35) provide much weaker constraints on the 
form of the Yukawa matrices than h = dzrj. In addition, we can obtain an arbitrarily small 
but hnite mass splitting AM^, independent of the Majorana mass term m/j, by choosing 
the Yukawa couplings such that 

Retrhmjj^h^ ~ Re tr ^ 77771 ^^^ 77 ^ + hm~^ri^ — . (4.36) 

This ability to tune the mass splitting of the light neutrinos independent of the magnitude 
of the Majorana mass term may have interesting implications in the light of the combined 
constraints provided by neutrino oscillation data and the current limits on lepton-flavour- 
violating and lepton-number-violating observables, including neutrinoless double-beta de¬ 
cay. 

A comprehensive phenomenological study of the aforementioned variations of this non- 
Hermitian Yukawa model in the context of current constraints from collider experiments 
and both astrophysical and cosmological observations (for recent reviews, see refs. [38, 39]) 
is beyond the scope of this article and will be presented elsewhere. 


5 Conclusions 

We have considered an extension of QED, whose non-Hermitian nature permits the sym¬ 
metry between the left- and right-chiral components of a Dirac fermion to be broken by the 
presence of an anti-Hermitian mass term. We have shown that the full gauge invariance of 
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this theory is restored when the Hermitian and anti-Hermitian masses are of equal magni¬ 
tude. Moreover, we have highlighted an intriguing possibility for explaining the smallness 
of the light neutrino masses and providing an additional source of CV violation through 
an analogous extension of the SM. Further phenomenological studies of this model and its 
variations are required in the context of the current low-energy neutrino data as well as 
both cosmological and astrophysical observations. 

Finally, we mention another direction of study, which deals with the dynamical genera¬ 
tion of the non-Hermitian mass term through nonperturbative quantum effects. Dynamical 
mass generation for neutrinos (with a vanishing bare mass) has been obtained in the con¬ 
text of Lorentz-symmetry violation [40-42], where the physical mass scale is provided by 
higher-order spatial derivatives. A nonperturbative mechanism could also be responsible 
for the non-Hermitian mass term in the present context, although the natural mass scale 
would be provided by the Higgs mechanism, instead of Lorentz-symmetry-violating opera¬ 
tors. In order to explore this avenue, one needs to derive a nonperturbative gap equation 
and study the possibility of a non-Hermitian mass term solution. Such a nontrivial solution 
could arise in a theory involving an axial coupling and is left for future work. 
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A One-loop corrections 


This appendix summarizes the technical details of the one-loop calculations described in 
Sec. 3. The elements A/j of a matrix A in the chiral field space are indexed by upper-case 
Roman indices I, J, K, M, ■ ■ ■ = L,R, where the LL element is in the top left. 


Passarino-Veltman parametrization. In d = 4 — 2e the two-point Passarino-Veltman 
form factors [33] are 


B. 


,;^;^,(p,mi,m2) = ^ 


1, kfj,, kp^kp 


— m\ + ie) ((p -I- kY — + ie) 

These may be related to the scalar form factors Hi, H 21 and B 22 via 


(AT) 


Bfj_{p,mi,m2) = Pf,Bi{p,mi,m2) , 

B^p{p,mi,m2) = p^iPuB2i{p,mi,m2) + r]fj,uB22{p-,rni,m2) 


whose divergent parts are 


Bo{p,mi,m2) D - , 


B2i{p,mi,m2) D — 


Bi{p,mi,m2) T - — , 

1 


B22{p,rni,m2) D {mi + m2 + ^ 


(A.2a) 

(A.2b) 


(A.3) 

(A.4) 
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In addition, we make use of the algebraic identities 

= ^[A{mi) - A{m 2 ) - {p'^ - mj - rnl)Bo{p,mi,rn 2 )] , (A.5a) 

p‘^B 2 i{p,mi,m 2 ) + dB 22 ip,rni,m 2 ) = A{m 2 ) + m\Bo{p,mi,m 2 ) ■, (A.5b) 

p^B 2 iip,mi,m 2 ) + B 22 {p,mi,m 2 ) = ^[A(m2) + {ml - ml - p‘^)Bi{p,rni,rn 2 )] , 


where A{m) is the tadpole form factor 


, fd<^k 1 

A{m) = / ^ -2 - 2 ”-— 

J iTT^ — m^ + le 


Lastly, for mi = m 2 we have the identity 


Bi{p,m,m) = - -Bo{p,m,m) . 


(A.5c) 

(A.6) 

(A.7) 

ivr^ (/j2 _ _|_p^^2 _ ^2 _|_ ( ) 

where the arguments of Co-^-,^u = Co-^-p,y{pi,p 2 -,mi,m 2 ,m'i) have been suppressed for 
notational brevity. We also define scalar form factors via 

= Pile'll + P 2 ^iCi 2 , (A.9a) 

= Pl^PmC '21 + P2 ^iP2uC22 + Pl{f,P2v)C22. + P^lvC2i • (A.9b) 

The only divergent form factor is (724, having the logarithmic divergence 


The three-point form factors are 

— / • 2 


1) fc/xj kfiki/ 


C 24 D ^ . 
4e 


(A.IO) 


Feynman rules. In the chiral basis, the Feynman rules of the model are [43-46] 
• To each photon line associate the factor (in the Feynman gauge) 




p^ + is 


(A.ll) 


To each chiral fermion line associate the factor 


iSij{p) = i 


.dijaj-p+ Mij 


(A.12) 


p2 — -|- is 

To avoid proliferation of sub- and superscripts, the spinor index assignment, denoted 
by the lower-case Gothic characters a and b, is understood as follows: 

• ^ab 


iSLL{p) = [iSLL{p\^ = ^2_M2 + ie 
iSrrIp) = [iSRR{p)]‘''^ = 

iSRiip) = [iSRL{p)r^ = 
iSlRip) = [iSLR{p)]a = 


ip ■ cr“^ 

p2 — -|- is ’ 

im-6°; 

__ b 

-|- is ’ 

im+6j’ 


p2 — -|- is ’ 


(A.13a) 

(A.13b) 

(A.13c) 

(A.13d) 
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with 


(A.14) 


= <>“ = <„ . 


-IJ. -u _ 

— a^’ 


• To each vertex associate a factor of — igijaj, where g = diag {g-,g+). 

• For any closed fermion loop include a factor of —1 and trace over the Lorentz indices. 

In the calculation of the one-loop corrections outlined below, we also make heavy use 
of the product and trace identities of the Pauli matrices, as listed in Appendix B of [46]. 

Fermion self-energy. The one-loop chiral fermion self-energies are given by 

/ d'^k 

ctk iSKNik + p) cr'^ iD^,y{k) , (A. 15) 

where we note that the couplings g± are dimensionful for d = 4 — 2e. The numerator is 


(2 — d) 6 kn ' k -\- dMxN- (A. 16) 

Rewriting in terms of the Passarino-Veltman form factors, we get 

^ij{p) = — d) dffiv CTic • pBi(p, M, 0) -|- dM kn Bo{p, M,0)] . (A.17) 

Hence, we obtain 


'Pll = 

d)cr ■pBi{p,M,0) , 

(A. 18a) 

'PRR = 

Xg^2 d)(T ■pBi{p,M,0) , 

(A.18b) 

P‘RL = 

107^2 ^ ^o{p,M,0) , 

(A.18c) 

'Plr = 

^ m+Bo{p,M,0) . 

(A.18d) 


The full fermion self-energy of the original Dirac field is obtained from the sum over 
the chiral indices I and J with correct weighting by projection operators. Specihcally, 


= PR'y^T,LLPL + PlJ^'PrrPr + PlPrlPl + Pr'PlrPr , 


(A.19) 


giving 




(A.20) 


Summing these contributions, we obtain the result in eq. (3.3). 
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Polarization tensor. The polarization tensor is given by 


/ Adh 

iSij{p + k) crXr iSNxik) ■ (A.21) 

Performing the trace over the Lorentz indices, we obtain the numerator 

2[(2k>"k'' +p^^k'''^ - - p^'^p ■ k)6ij 6 nk + Mjj M^k + ipuNK {p + k)f,kx] , 

(A.22) 

where is the Levi-Civita tensor. Here, we have dehned punk = 1, ifI = J = K = 

N = L, Punk = —1, iiI = J = N = K = R, and puNK = 0 otherwise. 

Rewriting in terms of the Passarino-Veltman form factors, we are left with 


n^"(p) = - 


9lK Qnj 


{ (pV - P^^P^) {B 21 + Bi)5u5nk 


- p^''[M‘^6uSnk - MuMnk)Bo 


+ ie 




PkP\{B21 + Bi) + Pk\p‘^B 22 PUNK^j ■ (A.23) 


When we sum over the chiral indices, the terms proportional to the Levi-Civita tensor 
cancel, and we obtain the result in eq. (3.8). 

Three-point vertex. The three-point vertices are given by 

/ d'^k 

iSuNik) iD^xik + p) iSpqik + p + q) aQ , 

(A.24) 

where p and q are the fermion momenta. The numerator is proportional to 


(2 - d) u'q a^k^ik + pd- q)p 5 kn ^pq + CTp Mkn Mpq 


+ 4 


k^ Skn Mpq + {k + p + q)^ Mkn Spq 


(A.25) 


such that the vertices can be written 


A^ - 


giK Qnp Pqj 
IGvr^ 


X I (2 — d) [{2 — d) ap C 2 A + o-q Up cr']^ Fi^p'^ 5 kn Spq + ap Mkn Mpq 

+ 4 (p^ Cii + q^ C 12 ) {6kn Mpq + Mkn ^pq) + 4 (p -|- q)^ Cq Mkn ^pq | ■ 

(A.26) 

Herein, we have dehned 

Fkp = PkPp iCii + C2i) + qK,qp {C22 ~\-Ci2)+PkQ.p iC2‘i + Ci\)-\-q^Pp {C2‘id-Ci2) ■ (A.27) 
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The three-point form factors are evaluated at = p, P 2 = Q, = w -3 = M and m 2 = 0. 
Hence, element by element, we find that 


AL = 

IGtt^ 

II 

cti 

9+ 

IGtt^ 

> 

II 

g+g- 


II 

g+g- 

dvr^ 


{2-d)\g.[{2-d)a^^C2A + dPa>^d^F^p) + g+d^^M^Co 
{2 - d)[g+(^{2 - d)a^^C24 + aPd^^a^F^p^ + g-a^M^Co 
m- {g+ + g-) {p^ Cn -b Cu) + g-{pP + qP)Co] , 

{g+ + g-)[pp Cn -b q^ Cn) + g+{l^ + q^)Co 


As in the case of the self-energy, we have 


AP = 


0 0 
1 0 , 


Kl + 


0 1 
0 0 


^RR + 


10' 

0 0 , 


Kl + 


0 0 
0 1 


^LR ■ 


(A.28a) 

(A.28b) 

(A.28c) 

(A.28d) 

(A.29) 


Summing over the contributions, we obtain the result for the total vertex in eq. (3.11). 
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